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Non-equilibrium diffusive systems are known to exhibit long-range correlations, which decay like the inverse 
1/L of the system size L in one dimension. Here, taking the example of the ABC model, we show that this 
size dependence becomes anomalous (the decay becomes a non-integer power of L) when the diffusive system 
' approaches a second-order phase transition. 

' This power-law decay as well as the L-dependence of the time-time correlations can be understood in terms 

of the dynamics of the amplitude of the first Fourier mode of the particle densities. This amplitude evolves 
^ ■ according to a Langevin equation in a quartic potential, 

' which was introduced in a previous work to explain the anomalous behavior of the cumulants of the current 

near this second-order phase transition. Here we also compute some of the cumulants away from the transition 
■ and show that they become singular as the transition is approached. 

^ ' Introduction 



X 



Ch ' Long range correlations are a well known property of non-equilibrium systems in their steady state UHZI]. For one dimensional 
(Jj , diffusive systems, which satisfy Fourier's law, these correlations (between pairs of points separated by a macroscopic distance) 
' scale like the inverse 1/i of the system size illi,|7|,|8i]. They can also be related to the non local nature of the the large deviation 
C/3 ', functional of the density profiles III ISD- 

^ ' These long range correlations have been calculated in a few cases, using in particular fluctuating hydrodynamics IH 0, S 
?3 ! It was observed that when a diffusive system approaches a second order phase transition, the factor in front of these 

^ i 1/L correlations becomes singular ifioll . The goal of the present work is to analyse these long range correlations at and in 
the neighborhood of a phase transition. To do so, we consider the ABC model, one of the simplest diffusive systems which 
undergoes a phase transition. 



Q ■ The main result obtained in the present work is that the long range correlations at the second-order phase transition of the 
O , ABC model have an anomalous dependence on the system size which can be understood by an effective theory ifTTIl for the 
' amplitude of the slow density mode which becomes unstable at the transition. 
I \ This effective theory, which was developed in a recent work ifTHl . led us to predict an anomalous size dependence of the 
■ cumulants of the particle current at the transition. Here, we also show that, away from the transition, these cumulants satisfy 
(NJ Fourier's law (i.e. are proportional to 1/ L), with prefactors which become singular at the transition. 

: The outline of the paper is as follows: in section |I] we first recall some known properties of the ABC model as well as our 
previous results! 1 1] on the cumulants of the current, which exhibit an anomalous Fourier's law at the second-order transition. In 
sections|II]and|IlIl we show how the effective theory for the slow density mode yields a L^^/^ size dependence of the correlation 
function in the critical regime. In section [TVl we calculate the cumulants of the particle current (all the cumulants in the flat 
phase and the first two cumulants in the modulated phase). They all decay like 1/L, with prefactors which become singular at 
the transition, and which match the expressions obtained directly in the critical regime in lUlll . 



I. SHORT REVIEW OF THE ABC MODEL 



' The ABC model is a one-dimensional lattice gas, where each site is occupied by one of three types of particles. A, B or C. 
Neighboring sites exchange particles at the rates 

AB ^ BA 
1 

BC ^ CB (1) 

CA ^ AC 

1 

with an asymmetry q < 1. This model has been studied on a closed and on an open interval (with particle reservoirs at each 
end) lfl2l - [l5ll as well as on a ring with periodic boundary conditions lfl6lj24ll : in this article, we consider the latter case by studying 
a ring of L sites. 
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For q = I, all configurations are equally likely in the steady state; on the other hand, for q < 1, the particles of the same 
species tend to gather iflql . When q scales as 

q = e-P'^ , 

the dynamics of the model becomes diffusive llzsl l26ll : for each species a G {A, B, C}, one can define a rescaled density profile 

Pa, 

Pro[site k is of type a] ~ pa{k/L, t/L^) = pa{x, r) . (2) 

Because the microscopic dynamics ^ conserves the numbers of particles, the Pa{x, r) are related to their associated currents 
ja{x, t) by the conservation laws 

dtPa{x,T) = -dx3a{x,T) . (3) 

By assuming that the underlying microscopic regions of the system are local equilibrium ifiol [Tsl IztIi . one can show that the 
currents ja{x, r) satisfy noisy, biased Pick's laws: 

ja = -dxPa + [iPaiPc - Pb) + -^r]a{x, t) (4) 

where h and c designate the previous and next species with respect to a, and where the ria{x, r) are Gaussian white noises such 
that 

{ria{x,T)ria'{x',T')) = (Taa'S{x ~ x')d{T - t') 

2pa{l-pa) ifa = a' 



with a, 



aa' 



-2paPa' Otherwise 



(note that these correlations imply riA + t/b + Vc =0 due to the fact that jA + in + jc is identically zero). At leading order in 
L, the noise in (|4]i can be neglected, so that (O leads to the evolution equations 11 1811 

drPa^ dlpa+ l3dxPa{Pb- Pc) ■ (5) 

For low /3, the flat density profiles Pa{x, r) = ra are a stable solution of these equations. They become linearly unstable lfl6[[l8[ 
I2TI1 when /? reaches /3* given by 

/3* = ^ with A-l-2^r2', 

indicating a second-order transition at /3 = /?* : above /3y , the steady-state profiles are modulated. It has been argued that these 
steady-state density profiles are time-independent ll2ll I23I1 . 

This stability analysis does not rule out the possibility of a first-order phase transition taking place at some /3t < /3* . A more 
detailed analysis of the neighborhood of shows that for A < 0, the transition should be first order ifisl . with 

^ = E'''-2E'a. (6) 



While no analytical expression for /3t is known in this case, it has been studied numerically in 12111 . 

In ifTHl . we considered the ABC model in the region A > where the transition is expected to be second order We found 
that, for a system of size L, there exists a critical regime |/3 — '-^ for which the dynamics of the density profiles is 

dominated by those of their first Fourier mode. Let 

1 T I / r f 1 if site k is of type A 

RA{t) = - V c-2-'=/^^fe t where Ak{t) = . (7) 

L I otherwise 

fe=i 

be the first Fourier mode of the density of species A. Then, near the transition, the first Fourier modes of the other densities, 
Rsit) and Rc{t), are related to Ra by 



RbH) = ^rc^^RAit) 
Rc{t) = '^--l+J^ RA{t) 



2r.4 (g) 
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and the evolution of Ra (t) can be described by a Langevin equation in a quartic potential in terms of the diffusive time t = t/Lp- 



dr 



with 



2A 



Ra 



(9) 



7 = 



/3* 



and with is a complex Gaussian white noise: 



{^,Aif)^,*Aif'))^24n^^-^^6{r-r') 



By rescahng Ra{t) by 



one can see that (|9|l becomes 



Ra{t) 



Ar\rBrc 



AL 



1/4 



fir), 



df 



(7-|/(r)P)/(f)+MT) 



where {tJ-{T)iJ.* {t')) = 6{t — r') and with 

I VSKrATBTc t 



A3/2 L5/2 



and 7 = \/L 



A3/2 



(10) 



(11) 



(12) 



Hence, in the critical regime |/3 — | ~ 1/ Vl, the amplitude of the first Fourier mode varies on a time scale f cx t/L^^"^, 
with an amplitude in 

In ifTltl, we showed that, due to these slow fluctuations of RA{t), the integrated particle current QA{t) of A particles dur- 
ing time t through a section in the system exhibits anomalous fluctuations at the transitions, reminiscent of those that can be 
numerically observed in momentum-conserving mechanical models ll28l : 

(13) 



(QaW) - \^PrA{rc - re) + JJJ^M 
{Q'Xit)) ^ A„, for n > 2 . 



We also argued that the coefficients An can be expressed in terms of n-point correlation functions of the solution /(f) of (fTTT i: 

1 



C„(7) = Jim -r / dfi..df„(|/(fi)../(f„)|2)^ . 



(14) 



Therefore in ifTlll we reduced the calculations of the cumulants of QA{t) in the critical regime to the study of the Langevin 
equation of a single particle evolving in a quartic potential. 

II. POWER-LAW RELAXATION OF THE FIRST FOURIER MODE 



In this section, we study the decay of the first Fourier mode in the critical regime when its initial amplitude is much larger than 
its steady-state values L^^^^). This is the case when one starts from a non steady-state initial condition: here we consider 
such a relaxation, starting from a fully segregated initial configuration of the type AA..AABB..BBCC..CC . 

According to [11], we expect all the higher Fourier modes of the densities to relax on the hydrodynamic time scale r = t/Li^. 
After this initial relaxation, the density profiles (|2]i should be dominated by a first Fourier mode with amplitude Ra [t] evolving 
according to (|9]l,(fT0]i: moreover, as long as Ra remains much larger than L^^/^, the noise term in (|9]l can be neglected and the 
evolution of Ra reduces to 



dr 



2A 



\Ra? 



Ra. 



(15) 
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When /3 = /3* (i.e. when 7 



0), Ra should thus decay as a power law: 

Ra{0) 



Ra{t) = 



1 



167rfA 



\Ra{0W 



Ra{0) a /77 
|i?A(0)| 47rV Ar 



(16) 



In Figure[T] we measured numerically the amplitude \Ra\ from its definition O for systems of 132 < L < 1500 particles for 
rA = rB ^ rc ^ 1/3 and /3 = /3*. The power-law decay (fTSTi should be valid in a rather limited range of time (L^ <C t <C L^^^), 
and is rather difficult to observe. Fortunately, the higher Fourier modes seem to relax fast enough for the power-law decay to 
occur already for t/L'^ ~ 10^^: our data for increasing sizes seems to converge to the power-law ( fTSI l in the whole range 

10-1 < t/L"^ < 2. 

Figure [U also shows a departure from this power law at a rescaled time t/L'^ increasing with L : in the next section, we show 
that this is due to the increasing effect of the noise term in (|9]l as Ra decreases. 

Power-law decay at the tricritical line 

The damping term of the evolution equation for Ra ( fTSl l that we obtained around (3 = P,, vanishes on the tricritical line A = 
(see eq. (|6]l). In this case, it is necessary to push the analysis of iflTll further in order to obtain an effective equation for Ra- One 
can show that, when A = 0, (flST l is replaced by 



dRA 

dr 



47r^ 7 



\Ra\' 
r2 A2 



Ra- 



(17) 



Along this tricritical line (7 = 0, A = 0), the t decay of Ra in the critical regime ( fTSl l should become in t ^Z^. 



in. CORRELATIONS IN THE STEADY STATE 



It has been shown that the ABC model exhibits long-range steady-state correlations, scaling as 1/L, in the flat phase /3 < /?* 
ifioll . In the modulated phase /?>/?*, these correlations should be of order 1 due to the modulation of the steady-state profiles. 

In this section, we show, using our effective dynamics (l9Tl,(fmi for the first mode Ra ©, that these correlations scale as 1/ VL 
in the critical regime; we also show that temporal correlations decay on the slow time scale f cx t/ L^^"^ at the transition. Finally, 
we briefly comment on the behavior on the tricritical line A = 0. 



10 



10- 



10- 



-1 _ 



10 



1 1 


1 


1 1 


L = 132 






L = 300 






L = 660 






L = 1500 






prediction 




1 1 


-1 




1 









FIG. 1: Relaxation of the first-mode amplitude measured using Q for systems of size 132 < L < 1500 started in a segregated initial 

configuration, for = t'b = re = 1/3 and P = 13,. Our analysis l |16l > of the first mode of the densities predicts a t^^^^ decay (see 1II6I1). 
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A. Spatial correlations in the critical regime 

As explained in iflTll . the density fluctuations in the critical regime are dominated by those of the first Fourier mode Ra- Thus 
one can calculate to leading order the steady-state correlations of the densities at the critical point: 



a a 

The Langevin equation (fTTT t can be expressed as a Fokker-Planck equation over the probabihty density of the rescaled first mode 
/(r), 

P(/,,/„r) = Pro[/(r):^/, + z/,]: 

dP 1 

— = div [(r^ - 7)Frl + -V^P (18) 
OT ■'4 

with f= {fx, fy) and 7 as defined in (fT2l l. From this equation, it is easy to see that /(f) is isotropically distributed in the steady 
state: 

Pq (r , r) cx exp \2^r'^ — r"'] . 

Therefore, {rI) ^ {RI^) ^ Q wA 



It is easy to compute 



J—CXD 



this leads to 



Pro[., = si]-Y,rlc. 2y^^J^Ci(7)cos27r^ . (20) 

Therefore, steady-state correlations should scale as 1/ \/L in the critical regime. In Figure HITAI we compare our prediction (|20b 
to numerical measurements of systems of 75 < ^ < 600 particles, for /? = and va = ''b = = 1/3- 



6 



I 



CO 

o 



10° 
10-1 
10-2 
10-3 
10-4 



: 1 


— L = 60 




L = 120 




L = 240 




L = 480 




L = 960 




- - prediction 


1 





0.1 



0.2 



FIG. 3: Steady-state time correlations Pro[sfc(t) — Sfc(O)] — X^r^ for systems of size 30 < L < 960, with ta ~ rs ~ rc = 1/3 and 
/3 = /?,, as functions of t/L^^'^ . Our effective, fluctuating evolution equation for the first Fourier mode ^ leads us to predict an exponential 
decay ( 121b . 



Remark: The 7 — >^ —00 and 7 — > +00 limits of our expression for the equal-time correlations ( l20l i can both be 
checked from known results. 

In iToll . exact expressions for the equal-time density correlations have been calculated in the flat phase ( fioll . eq. (23)): they 
diverge as /3 — ^ leading to 

Pro[., ^s,]-^ ^ ^—^ cos 2.-^ , 

which is compatible with 6*1(7) — — l/22_for 7 — )• —00 in ( fT9] l. 

On the other hand, Equation (39) of iflSll establishes that the modulated steady-state profiles in the disordered phase are, for 
/3 -> of the form 



with X = k/L and ip an arbitrary phase. This leads to 



1 v-2 A2/?-/3, ^ k-l 
Pro[sfe = si] - 2^r„ ~ — cos27r— ^ 



which is compatible with Ci (7) ~ 7 for 7 ^ cx) in ( fT9] l and ( l20l i. 



B. Decay of the steady-state temporal correlations 

In the previous section, we have seen that, on a time scale ^ t <C L^^^, the first Fourier mode of the densities ^ should 
decay as a power law when the system relaxes from a non steady state initial configuration. Here we consider the decay of 
the steady-state density correlations {pa{x, 0)pa{x, t))^: the noisy evolution equation (|9]l implies that they should decay on the 
"slow" time scale t oc L^/^. They can be expressed in terms of the rescaled first mode /(f) as 

Pro[s,(0 = Sfe(O)] - ^ (X (/(0)/*(f ) + cc), 

with f cx r/VX the rescaled time (fT2l) . Let be the operator of the Fokker-PIanck equation dTSI) over the probability density 

of/(f),P(/,,/„f): 
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For large r, correlations such as (/(O)/* (r)) ^ should decay as e/^'''^ , where A^^ is the second largest eigenvalue of (the largest 
being 0). Thus, we expect time correlations to decay on the time scale t cx L'^/^ as 

t 



Pro[sfc(0) = Sk{t)] - ^ cx exp 



(21) 



with a — —Stt^ ^^^^^^^ \^ . While there is no known analytical expression for A^, one can determine it numerically by 
approximating the operator over a finite subspace of L^(R^) of growing dimension. Because the steady-state (and thus the 
0-eigenvector of H^) is Po{fx, fy) oc e^^^'~^ , we consider the subspace 

where Q{fx, fy) is a polynomial of degree less than some finite N in (fx, fy). We need to define a scalar product (-I-) for which 
Hj is Hermitian: to do so, we choose 

{Pq\Pb) =11 dxdye^"-^^''" Pq{x, y)PR{x, y) (22) 
dxdye^'^''^~''''Q{x, y)R{x, y) , 



for which the matrix elements of read 



1 



{Pq\H^\Pr) ^-ijj dxdye^^'^"-'\^Q).{^R) . (23) 

For 7 = 0, we constructed an orthonormal basis of the subspace of the Pq with respect to the scalar product (l22l l for 2 < < 8; 
then, we computed the matrix elements of Hq ( |23] l in this basis. The second largest eigenvalue of these successive hermitian 
matrices appears to converge to 

Ao ~ -0.83. 

In Figure HITBI we compare this prediction with the results of simulations in the steady-state for 60 < L < 960, = rs = 
rc = 1/3 and/? = 

The exponential decay (|2T]) also seems to predict the cut-off of the power law decay of Figure[T]observed in systems relaxing 
from an arbitrary initial condition to the stationary state. 



C. Correlations on the tricritical line 



As in the deterministic case the cubic term of the fluctuating evolution equation for the first Fourier mode (|9]l vanishes on the 
tricritical line A = 0. The fluctuating correction of the deterministic tri-critical evolution equation STl\ is the same as those of 
® : 



dRA 



\Ra\' 



= 47r^ 7 - „ ^ 



Ra 



Vl 



This leads to the rescaling 



with 



= 9{r) with — = (7 - \g\ )g + ^(r) , 



47r2 



,2/3. 



L8/3 



and 7 



Hence, the tricritical regime occurs for |/3 — oc L~^^^ : it is characterized by fluctuations of the first Fourier mode of 
amplitude on a time scale t cx L^^^. 

In its steady state, we expect spatial correlations to scale as 

^ 1 k — I 

Pro[sfc(t) = Slit)] - 2]^a °^ ' 

while time correlations should decay exponentially on the time scale t cx L^^^, as 



1 



Pro[sA;(G) = Skit)] - cx j-^ exp 



t 
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IV. FLUCTUATIONS OF THE CURRENT IN THE ABC MODEL 

In this section, we study the cumulants of the integrated current of particles of type A, QA{t), as /3 — > 

In contrast to ifTHl . where we studied the critical regime |/3 — /3* | ^ 1/ \/L, here we compute the cumulants to leading order 
in L for a fixed /3 7^ /3* : the expressions we obtain diverge at /?* . For /3 — ?> /3~ , we find (|36] | 

= * + andMW).-, 2= ,24) 




while we find that, for -> /3+ (l34] i. 



2 



We also show that these expressions are compatible with the 7 — >^ ±00 limits of the expressions of the cumulants (iJTl l.dSSll in 
terms of the functions C„(7) ( fT4l i derived in the critical regime in ifTTl] . 

A. First cumulant of the current 

We start with the assumption that the steady state density profiles Pa{x), which are the long-time limits of the leading-order, 
deterministic hydrodynamic equations (|5]l, are time-independent, both in the flat and in the modulated phase. The associated 
particle currents ja are thus homogeneous, ja{x, r) = Ja, and can be computed by integrating ^ in space, yielding 

Ja= I I3pa{x){pc{x) ~ Pb{x))dx . (25) 
Jo 

Then, the average integrated current {QA{t)) through any position in the system will behave in the long-time limit like {QA{t)) — 



;|;Ja, so that 

(Q 

t LJo 

For /3 < we thus obtain {QA{t)) — j^PrA{rc — fs)- For /3 > B^, the analytic expression of the Pa{x) is rather 
complicated 112 ill . However, their limit as /3 — > /3+ takes a simple form (see llST . or the 7 — > oo Hmit of dill): 

Pa{x) ~ r, + (i?,c2''^" + cc.) (26) 



with \Ra\ = ^\/'rA{P^li*)/ 2 A/3* and with Rb, Rc related to Ra by (O. This leads to an analytical expression for {QA{t)) 
when /3 /3+ : 



{QA{t)) _ rc-re 
t PiP, L 
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B. Second cumulant of the current 



In order to compute fluctuations of the current QA{t)^ the noise terms of the biaised Pick's law (HJi have to be taken into ac- 
count. These fluctuati ng hydr odynamics can be reformulated as a large deviation principle known as the macroscopic fluctuation 
theory (MFT)|[i3|2l&l3ll], which gives the probability of observing a given evolution Pa{x, r) of the density profiles : 



Pro[pa(a;, r)] oc cxp 



-L II dxdT^ij ~ q).<T ^{j-q) 



(27) 



where j = {JaiJe), q = {qA,qB) (with qa = —dxPa + PPa{Pc — Pb) the deterministic part of Q) and with 

Pa(I-Pa) -PaPb 
~PaPb - Pb) 



f ctaa 


f^AB 1 


-A 


\ cfba 


<7BB J 
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The generating function log ^e-^Q'4(*)^ can then be calculated as the solution of a variational problem : 



log(e 



,AQ.4(t) 



maxL / dr I dx 

P'^da Jq Jq 



(28) 



For A = 0, the optimum in the above is the solution to the deterministic equations (|5]l, which satisfy by definition j ~ q '■ these 
are pa{x, t) = (for /? < P^) or pa{x, r) = pa{x) (for /3 > and ja{x, t) = Ja defined by 

For (3 < j3^,, the optimal profiles are still flat even for A 7^ (as we will see below, they are a stable minimum); 



p{x,t) = r 

j{x,T) = «7 + (T Q j = J(A) 

with p ~ {pA, pb), r = {vA, tb), and q and cr taking their constant values for p — r. This leads to 
and to 



A 



~ 4 [/3ArA(rc - tb) + A2m(1 - ta)] = ^F^M) 



(29) 



— ryi(l — r^) and — - = I ^ I for n > 3 . 



(30) 



(31) 



On the other hand, the optimal profiles in ( |28] | vary with A for with non-trivial optimization equations. We will 

therefore restrict ourselves to the calculation of {Q\{t)), for which it is sufficient to compute ( |28] ) to second order in A. 

In order to do so, we consider profiles close to the A = optimum, pa = Pa{x) and ja — Ja- Because arbitrary translations 
of these profiles are also optimal, we need to consider profiles moving at a small velocity v : this amounts to 



p{x, t) = p{x — vt) + p{x — vt) 
j{x,t) = J + K — vp{x — vt) 

with p = {pa, Pb) and yu(x), v,K Because qa = —dxPa + Ppa{pc — Pb), this leads to the following expression for q : 



q = q- fi - Mp with M = /3 
q — q{p) and p' = dxP- Then, the right-hand side of ( |28] ) becomes 



1 - 2pc 2pA 
-2pB 2pc-l 



S = max — / dx \ \{Ja + Ka — vpa) 



^{K -vp + p' + Mp).a~^{K -vp + p' + Mp) 



with (T = <t[p). The optimization equations over K, v and p can be written as 

A 



dxF 







dxp.F = Xta 
F' = M^F 

with F — (j^^{K — vp + p' + Mp). The right-hand side of ( |28] | then reads 



(32) 



5* = max — I dx 

K,v,F{x) L 



\{Ja + Ka~ vpa) - -F.aF 
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It can be expressed completely in terms of the optimal F{x) solution of (|32] | by using 



j dxF.aF = J 



F.{K - vp + ^i' + M^i) 
K. / dxF~v / dxp.F+ / dxfi.{M^F - F') 



= X{Ka - vta) 
so that the generating function is given to second order in A by 



log(e^Q-W)~i J dx 







j dx 


XJa + \f.ctF 



(33) 



By determining the solution F{x) of (|32] | and integrating (l33T l numerically, {Q\{t)) ^ can be predicted to leading order in L for 
/3 > /3* . This prediction diverges as /? /3+ : because p{x) is known analytically in this limit ( |26] ). we were able to calculate 
F{x) exactly, 



F{x) 



AnX^rs ~ rc) 



A3/2(/3-/?,) \-Ra 
leading to a simple expression for in this limit: 

t 127rr^rBrc(rB - rp)^ 



cc. + 0(1) 



(Qlit)) 



AA(/3 - 13., 



(34) 



C. Higher-order corrections in the flat phase 



As shown above, the macroscopic fluctuation theory ( |28] l predicts, for large L, Gaussian fluctuations of Qa [t) in the flat phase 
/3 < /3*, with a non-singular variance dsTl ) as /3 — )■ /3~. 

Here, we calculate the next 1/ L corrections to the generating function ( l30l l. by generalizing to the ABC model the approach 
followed in lf32ll in the case of a single conserved quantity. These corrections, of order l/i^, can be obtained by considering the 
large-deviation principle (|27] | as a functional integral : 



« j D[pa]D[ja]cyi^ L jj dxdr yXjA-\{i - q)-(y ^(j-q) 



(35) 



where the integral takes place over all profiles {paiia) compatible with the conservation law dxja + drPa = 0. Fluctuations 
around the optimum profile ( |29] l then give corrections to the saddle-point expression ( |30] l. We now consider such fluctuations, 
expressing them in terms of their Fourier modes : 



p{x, t) = r + 5p{x, t) 
j{x,T) = J(A) +6j{x,T) 



with 



6pix, r) = ^ fc [ak^e'^"^-'^'^ + cc. 

k.uj 

6j{x, r) = ^ w [afc^e'(^^-'^*) + cc. 



with a^uj ~ {(^kuj ' "-"^iw^ ) amphtude of the fluctuations of wave number k and pulsation uj, which take discrete values : 

27rm 

k = 27rn with n e N and u ~ , „ with m ^"L . 



Expanding ( |35] | to second order in the aki^, we obtain 



Y^dakujdal^ 


cxp 


k.LU 





L 



Ffiat(A) 



k.LU 
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with i^flat (A) the dominant-order generating function ( |30] | and 

(j)A + (t)c (t>c - 6i/3fc3 
0c + 6i/3/c3 (j)B + (t)c 



M, 



, with 



Integrating the Gaussian variables akuj gives the corrections 

log (e^«- ) ^ I Fflat (A) - 5] log (A, c.) + C 



with Pfe(A, oj) = + 4>b4'c + 4'c4'A — 36/3^fc^ and C an additive constant fixed by the condition log ^e'*''^') = for A = 0. 

In the t oo limit, the sum over w can be replaced by an integral : 



|Pflat(A) - ^ 



f dLj\ogPk{\u:)+C. 

1.. 'J 



Because only the first mode of the fluctuations, k = 27r, becomes unstable as /3 — >■ we expect the divergence in the 
cumulants at the transition to only affect this k ~ 2tt term. Thus we take the following limit in P2tt (A, w) : A ~ (to determine 
the cumulants), w ~ (we expect slow fluctuations to be responsible for the divergence), and /3 ~ Then P27r(A, uj) takes the 
simplified expression 



P2,(A,c<.)~64^^ 



w2+4^2^(/3_/3^)2 



rArsrc 

The integral of log Q{X,ll!) is apparently divergent : however, since 



+ 12/33(rc-rB)A 



= Q(A,Lj). 



dajlog(5(A,a;) 



4n log n 



12k 

- 47r^\/ A(/3 - (3ky + -^rArBrcirs - rc)X , 



its divergent part is canceled out by adjusting C so that log ^e'^'^' ) = for A = 0. Hence, the part of the generating function 
which becomes singular as /3 /3» reads 



log(e^<3^(*)).lp,a,(A) + l^ 



247r rATBrc 



leading to a divergence of the n-th cumulant of Qa (t) scaling as ( l24b . with 



(rs - rc)X - 1 



Br, 



-r(n- 1/2) 



2ATTrArBrc{rc ~ rs) 

A572 



(36) 



D. Anomalous fluctuations in the critical regime 



In ifTHl . we found that the fluctuations of the first Fourier mode on the slow time scale r cx t/ L^l"^ ( fT2] ) lead to anomalous 
fluctuations of the integrated current of QA{t) in the critical regime (fT3] l. More precisely, we derived 



t 2t 



t Stt^VSA 



L5/2- 



A3/2 



A^/^{rb - r, 



SAvanrc 



A 



Ci(7) 



27rA 



(37) 
(38) 



with 7 and Cni'j) as defined in (fT2l) and (fT4l i. Because /(f) evolves in a quartic potential (fTTl i. only Ci(7) can be easily 
calculated ( fT9] ). However, (|37] ) and ( |38] | predict the dependence of the cumulants in {rA,rB,rc, /3) in terms of the unique 
parameter 7 ( fT2b : one can easily check that ( l36b and ( |34] | are consistent with this dependence, for 
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V. CONCLUSION 

In this paper we have shown that the long-range correlations ( l20l i of the ABC model near the second-order phase transition 
decay like the L^^/^ power of the system size L. In the entire critical regime, these correlations UtI l23ll28ll can be understood 
from the evolution (11] of the amplitude of the first Fourier mode given by the Langevin equation of a particle in a quartic 
potential ( fTTT i. 

We have also computed the cumulants of the current of particles jTM away from the transition, showing that the become 
singular at the transition in a way which matches the results of our previous work ifTlll where these cumulants were computed in 
the critical regime. 

It would be interesting to see whether other diffusive systems, at a phase transition, display correlation functions and current 
fluctuations with behaviors similar to those we discovered here for the ABC model. 

Deterministic one-dimensional systems, in particular those which conserve momentum, are known to exhibit an anomalous 
Fourier's law ll33ll . with cumulants of the current ll28ll and correlations ll34l - l36ll scaling as a non-integer power of the system 
size. Although these systems are much more difficult to study than the ABC model (for which one only needs to follow the 
dynamics of a single mode), it would be interesting to see whether the approximations that have been used so far, such as the 
mode-coupling approach 1137143911 . could predict the power-law dependence of these current and density fluctuations. 

Acknowledgments. BD acknowledges the support of the French Ministry of Education through the ANR 2010 BLAN 0108 01 
grant. 
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